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THE  BOX  METHOD  FOR  LINEAR  PROGRAMMING: 
PART  II  -  TREATMENT  OF  PROBLEMS  IN  STANDARD 
FORM  WITH  EXPLICITLY  BOUNDED  VARIABLES 


by  Karel  ZIKAN  and  Richard  W.  COTTLE 


1.  Introduction. 

The  Box  Method  is  a  new  interior-point  algorithm  for  linear  programming  that  deals  with 
linear  programs  whose  constraints  are  all  linear  inequalities.  [See  Zikan  and  Cottle  (1987, 
Algorithm  I)  .]  Hence,  from  the  standpoint  of  (primal)  problems  presented  in  what  Dantzig 
(1963,  p.  86)  calls  “standard  form”  (i.e.,  equality  constraints,  nonnegative  variables),  the 
Box  Method  operates  on  the  dual  problem;  in  this  sense  it  is  a  dual  method. 

Since  linear  programming  problems  with  explicit  upper  (and  lower)  bound  constraints 
are  very  common,  a  practical  linear  programming  algorithm  must  be  able  to  treat  them  as 
easily  as  the  Simplex  Method  does.  This  report  aims  to  accomplish  for  the  Box  Method 
what  Dantzig  (1955)  did  for  the  Simplex  Method:  to  produce  an  efficient  variant  of  the 
algorithm  for  handling  problems  with  (explicit)  upper  bounds  on  the  variables.  The  ex¬ 
plicitly  bounded  variables  are  assumed  to  occur  in  a  primal  problem  having  standard  form, 
although  the  algorithm  will  solve  its  dual.  Like  the  Simplex  Method,  the  variant  presented 
here  solves  the  problem  almost  as  efficiently  as  if  the  upper  bounds  were  not  present.  The 
increases  in  problem  size  and  computational  effort  are  very  small. 

Specifically,  it  will  be  shown  in  this  paper  how  a  crucial  aspect  of  the  Box  Method 
can  take  advantage  of  the  special  structure  of  the  dual  of  a  linear  program  in  standard 
form  with  explicit  upper  bounds  on  the  variables.  The  key  step  in  the  Box  Method  is 
the  determination  of  a  minimum-weight  basis  corresponding  to  an  interior  feasible  iterate. 
(This  is  also  called  the  combinatorial  or  matroidal  problem.)  Minimum- weight  bases  are 
used  in  constructing  the  “boxes”  over  which  the  algorithm's  direction-finding  subproblems 
are  solved.  Since  the  other  algorithmic  effects  are  insignificant,  only  matters  related  to  the 
minimum-weight  basis  will  be  discussed.  An  important  consequence  of  the  results  presented 
here  is  that  any  efficient  scheme  pertaining  to  (primal)  problems  without  explicit  upper 
bounds  on  the  variables  can  be  extended  to  the  bounded-variable  case.  This  point  will  be 
illustrated  in  the  sequel  to  this  report  where  the  Box  Method  will  be  specialized  to  the 
important  class  of  (linear)  minimum  cost  flow  problems. 


2.  Block  formulation. 
Consider  the  linear  program 


minimize  b (1.1) 
subject  to  A^y  =  c  (1.2) 

0  <  y  <  u  (1.3) 

where  ATis  annxm  matrix  with  linearly  independent  rows.  Accordingly,  rank  Ar  =  n  <  m. 
After  slack  variables  are  inserted  in  (1),  the  constraint  matrix  becomes 


the  objective  function  vector  and  right-hand  side  vector  become,  respectively, 

F=(»T,0’)  «,d  *=(o)- 

It  is  immediate  that  AT  is  an  (n  +  m)  x  2m  matrix  of  full  rank.  The  purpose  of  an  upper- 
bounding  technique  is  to  avoid  using  such  a  potentially  large  matrix  and  to  work  instead 
with  AT  itself. 

Now  assume  that  the  columns  of  the  matrix  A7  are  “weighted”  as  in  the  Box  Method.1 
This  means  that  to  each  column  of  AT  is  associated  a  positive  real  number,  its  weight. 
How  the  weights  are  assigned  is  not  important  at  the  moment.  Given  the  weight  vector, 
w,  the  first  problem  is  to  find  a  minimum-weight  basis  consisting  of  columns  of  the  matrix 
AT  and  then  an  appropriate  factorization.  A  minimum-weight  basis  is  a  maximal  set  of 
linearly  independent  columns  for  which  the  sum  of  the  corresponding  weights  is  minimum. 
It  will  be  shown  that  solving  this  matraidal  problem  MP(A1 ',  w)  can  be  accomplished  with 
very  little  effort  beyond  what  is  needed  to  do  the  analogous  job  for  AT.  Furthermore,  the 
structure  of  a  minimum-weight  basis  facilitates  solving  the  two  equations  for  which  the  Bax 
Method  requires  solutions. 

3.  Theory. 

The  propositions  and  corollaries  of  this  section  provide  the  theoretical  underpinnings  of  the 
algorithmic  developments  that  follow.  The  results  proceed  from  the  natural  (and  obvious) 
pairing  between  the  columns  of  the  (n  +  m)x  2m  matrix  AT,  namely  column  t  corresponds 
to  column  m  +  *.  Therefore,  let 

i  « — ►  i'  =  m  +  i  *  =  l,...,m.  (2) 

‘The  Box  Method,  as  stated,  would  deal  with  A  (rather  than  A7)  and  hence  the  row*  would  be  weighted. 


The  first  result  is  for  matrices  of  more  general  type  than  Xs. 

Proposition  1.  Let  M  be  akxt  matrix  having  no  column  of  zeros.  Assume  B  and  D  are 
k  x  k  invertible  and  lx  l  invertible  diagonal  matrices,  respectively.  Given  a  fixed  £- vector 
of  weights,  tv,  the  matroidal  problems  associated  with  the  matrices  M  and  M  =  BMD  are 
solved  by  the  same  set  of  indices. 

Proof.  The  matrices  Af  and  BMD  have  exactly  the  same  sets  of  linearly  independent 
columns.  That  is,  if  E  denotes  a  nonempty  index  set  drawn  from  {1, . . . ,  f},  then  M.e  has 
linearly  independent  columns  if  and  only  if  M.e  has  linearly  independent  columns.  Since 
the  two  matroidal  problems  axe  defined  in  terms  of  the  same  weight  vector,  tv,  it  follows 
that  the  greedy  algorithm  will  make  the  same  column  choices  in  each  problem,  thereby 
producing  the  same  index  set  as  solution.  ■ 

Corollary  1.  Assume  that  the  matrix  Ar  and  the  vector  of  weights,  tv,  are  fixed.  Then  the 
same  set  of  basic  indices  solves  the  matroidal  problem  for  every  matrix  having  the  block 


where  B  is  any  mxm  invertible  matrix  and  D  is  an  invertible  diagonal  matrix  of  the  same 
order. 


Proof.  In  view  of  Proposition  1,  it  suffices  to  note  that  a  matrix  of  the  form  (3)  can  be 
brought  to  the  "canonical  form”  (2)  through  multiplication  on  the  left  by 


and  on  the  right  by 


(o  B°") 
(o  D")' 


In  the  case  of  (4),  the  identity  matrix  is  of  order  n  whereas  in  (5)  it  is  of  order  m.  ■ 

Remark.  Assumption  (A3)  of  the  Bax  Method  ordinarily  calls  for  some  (row)  scaling  in 
the  A.  Accordingly,  one  is  required  to  solve  matroidal  problems  in  which  the  matrices  have 
the  block  form 


/  AtD  0  ) 

V  D  I)' 


but  such  matrices  are  easily  reduced  to  the  form  (2). 
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Corollary  2.  Assume  the  matrix  AT  and  the  vector  of  weights ,  to,  are  fixed.  Then  the 
same  set  of  basic  indices  solves  all  matroidal  problems  for  matrices  of  the  form 


where  D\  and  D?  are  nonsingular  diagonal  matrices  of  order  m. 

Proof.  Obvious.  ■ 

Remark.  A  special  case  of  Corollary  2  is  that  in  which  D\  and  are  (possibly  distinct) 
"signed  identity  matrices”  (i.e.,  their  diagonal  elements  are  ±1.) 

Tie-breaking  assumption.  In  order  to  avoid  clumsiness  in  the  wording  of  the  following 
results,  it  will  be  assumed  that  ties  are  broken  according  to  an  arbitrary — but  fixed — rule. 
This  being  settled,  no  special  consideration  need  be  given  to  ties. 

Proposition  2.  fYom  each  pair  of  indices  (t,  *'),  at  least  one  index,  namely  the  one  with 
the  smaller  associated  weight,  must  be  selected  when  the  greedy  algorithm  is  applied  to 
the  matroidal  problem  associated  with  AT. 

Proof.  Given  the  special  structure  of  the  matrix,  the  column  »'  can  not  be  linearly  depen¬ 
dent  on  previously  selected  columns  unless  its  "mate”,  »,  has  already  been  selected.  The 
same  holds  for  column  i.  Since  the  greedy  algorithm  considers  the  columns  according  to 
increasing  order  of  their  weights,  the  assertion  readily  follows.  ■ 

Proposition  3.  All  bases,  i.e.,  maximal  linearly  independent  sets  of  columns,  in  the  matrix 
AT  contain  exactly  n  (i,  i')  pairs  and  m  —  n  single  elements  from  the  remaining  ( » , «')  pairs. 

Proof.  The  rank  of  AT  is  n  +  m,  hence  this  is  number  of  columns  in  any  maximal  linearly 
independent  set.  Let  E  denote  the  set  of  indices  of  such  a  set  of  columns.  There  are  m  pairs 
of  indices  (i,  i'),  and  E  must  contain  at  least  one  representative  of  each  pair.  (Otherwise, 
the  corresponding  matrix  would  contain  a  row  of  zeros  and  hence  would  not  consist  of 
linearly  independent  columns.)  Assume  these  m  indices  have  been  identified.  No  matter 
how  the  other  n  indices  are  chosen,  their  mates  will  necessarily  belong  to  E.  After  this 
pairing  is  established,  the  there  must  remain  m  —  n  among  the  original  m  elements  for 
which  E  contains  no  mate;  these  are  the  single  elements.  ■ 

The  idea  is  illustrated  by  the  matrix 

f  Aj  Aj  0  0  \ 

(/>,),  0  (D,),  0 

\  0  (D\)2  0  (f?j)a  / 
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in  which  Aj  denotes  a  nonsingular  n  x  n  matrix;  for  *  =  1,2,  (D<) j  and  (Dfo  denote 
nonsingular  diagonal  matrices  of  orders  n  and  m  —  n,  respectively.  In  the  present  case, 
the  pairs  of  indices  would  be  (1,1'),...,  (n,n');  this  corresponds  to  the  fact  using  Aj  to 
cover  the  first  n  rows  makes  it  necessary  for  (Da)i  to  cover  the  next  n  rows.  The  remaining 
rows  must  be  covered  by  columns  drawn  in  a  “complementary”  way — that  is,  exactly  one 
of  *,t' — from  the  matrices  (1)1)3  and  (1)3)3. 

Definition.  Given  w  €  i?2"*,  let  w  €  A3*”  be  defined  by  the  relations 

Wi  =  max  {u\,  u>m+,},  wm+i  =  min  {u>,,  (8) 

The  following  result  pertains  to  the  selection  of  a  basis  in  accordance  with  the  associated 
matroidal  problem. 

Proposition  4.  For  a  fixed  matrix  AT  the  solution  E(w)  of  the  matroidal  problem 
MP(AT,  tv)  associated  with  the  weight-vector  w  and  the  solution  E(w)  of  the  matroidal 
problem  MP(AV,  w)  associated  with  with  weight-vector  w  contain  the  same  set  of  pairs  of 
indices. 

Proof.  Consider  the  problem  MP(Ar,  w)  associated  with  the  matrix  Ar.  Assume  w  ^  w, 
for  otherwise  there  is  nothing  to  prove.  Define  the  set 

U  =  {i  €  {l,...,m}  :  wt  >  to,} . 

Let  t  €  U  be  arbitrary,  and  perform  the  following  sequence  of  operations: 

1.  Multiply  AT  on  the  left  by  the  elementary  matrix  that  corresponds  to  a  Gauss- Jordan 

pivot  on  the  element  in  row  n  +  i  and  column  »,  that  is,  (Cf.  Proposition  1: 

left  multiplication  by  nonsingular  B.) 

2.  Change  the  sign  of  column  i'  =  m  +  i.  (Cf.  Proposition  1:  right  multiplication  by  a 
nonsingular  diagonal  matrix.) 

3.  Interchange  columns  t  and  t'  and  also  interchange  the  corresponding  pair  of  weights. 
(Otherwise,  the  matroidal  problems  are  not  the  same.  Note  that  a — nondiagonal — 
permutation  matrix  is  being  applied  on  the  right!) 

4.  Bring  the  matrix  (which  has  the  form  (7)  with  signed  identity  matrices)  to  the  form 
(2).  (Cf.  Corollary  2.) 

Now  the  structure  of  the  matrix  is  the  same  as  in  the  beginning  except  for  the  roles  of 
the  indices  i  and «'  which  have  been  interchanged.  Repeat  the  same  procedure  for  all  t  €  U. 


When  the  process  is  complete,  the  original  matroidal  problem  MP(AT,  to)  will  have  been 
transformed  into  the  canonical  format  of  the  matroidal  problem  MP(AT,  to).  The  difference 
is  that  for  each  t  €  U ,  the  column  in  the  i  position  now  bears  the  name  t-/  and  vice  versa. 
Therefore,  for  each  t  6  Ut  i  6  E(w)  if  and  only  if  i'  6  E(w),  and  %'  €  E(w)  if  and  only  if 
*  €  E(w).  Likewise,  for  each  «  ^  U,  i  €  E(w)  if  and  only  if  t  €  E(w\  and  *'  e  E(w)  if  and 
only  if  V  €  E(w).  Thus,  E(w)  and  E( to)  contain  the  same  index  pairs.  ■ 

Corollary  3.  To  solve  the  matroidal  problem  MP(AT,w),  it  suffices  to  solve  MP(Ar,  to) 
where  to  denotes  the  first  m  components  of  the  vector  to. 

Proof.  Obvious.  ■ 

Example.  Suppose 

^T=  (  1  j!)  1111(110  =  (5>8, 6, 9, 4, 7), 
so  that  m  =  3  and  n  =  2.  Thus, 

to  =  (9, 8, 7, 5, 4, 6)  and  U  =  {1,3}. 

To  solve  the  matroidal  problem  MP(AT,  to)  one  must  select  5  linearly  independent  columns 
from  the  matrix  ATfor  which  the  sum  of  the  corresponding  weights  is  as  small  as  possible. 
The  data  can  be  arranged  in  tabular  form  as  follows: 


1 

2 

3 

V 

2! 

3' 

< —  Indices 

5 

8 

6 

9 

4 

7 

< —  Weight  vector  to 

5 

2 

3 

0 

0 

0 

1 

2 

4 

0 

0 

0 

1 

0 

0 

1 

0 

0 

<—  AJ 

0 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 

The  transformation  process  described  in  Proposition  4  ends  with  the  table 


1' 

2 

3* 

1 

2' 

3 

< — 

9 

8 

7 

5 

4 

6 

< — 

5 

2 

3 

0 

0 

0 

1 

2 

4 

0 

0 

0 

1 

0 

0 

1 

0 

0 

< — 

0 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 
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Indices 

Weight  vector  to 

A 1 


1 

« 

I 


The  solution  to  MP(AT,w)  is 

E(w)  «  {2,3',1,2',3}. 

By  contrast,  the  matroidal  problem  MP(AT,  w)  is  expressed  by  the  table 

< —  Indices 
< —  Weight  vector  w 

«—  AT 

which  only  differs  from  the  preceding  table  in  the  locations  of  the  primes  on  the  indices 
(i.e.,  the  columns  labels).  The  solution  of  MP(Ar,w)  is  obviously 

E(w)  —  {2,3,l',2/,3'}. 

Proposition  2  asserts  that  for  every  i  =  1, . . .  ,m,  the  index  of  the  smaller  weight 
belongs  to  E(w).  The  definition  of  tv  implies  that  *'  €  E(tb)  for  every  t,  and  then  focuses 
attention  on  solving  MP(Ar,  w)  where  vbi  =  tp,-  for  i  =  1, . . . ,  m. 

4.  Algorithm. 

As  seen  above,  the  solution  set  E(w)  contains  n  pairs  of  indices  (t,  i')  and  m  —  n  other 
elements,  one  from  each  remaining  (t,  t')  pair,  namely  the  one  corresponding  to  the  smaller 
weight.  Accordingly,  the  set  E(w)  has  a  natural  partitioning  into  the  disjoint  union  of  four 
subsets:  D,  D',  5,  and  S'.  The  following  relationships  hold. 

i)  D  and  5  are  subsets  of  {1, . . .  ,m); 

ii)  D'  and  S'  are  subsets  of  {1', . . .  ,m'}; 

iii)  *  €  D  if  and  only  if  *'  €  D*\ 

iv)  D  and  U  each  contain  n  elements; 

v)  5  and  S'  may  each  be  empty,  but  their  union  contains  exactly  m  —  n  elements. 


1  2  3  V  2!  2? 

9  8  7  5  4  6 

5  2  3  0  0  0 

1  2  4  0  0  0 

10  0  10  0 

0  10  0  10 
0  0  1  0  0  1 


procedure  upper  bounds  min- weight  basis 
begin 

S  < —  0;  (comment:  S  is  the  set  of  single  t’s  selected) 

S' «—  0;  (comment:  S*  is  the  set  of  single  t'’s  selected) 

D  *—  0;  (comment:  D  is  the  set  of  t’s  selected  for  pairs) 

XK  •* —  0;  (comment:  D'  is  the  set  of  t'’s  selected  for  pairs) 

for  t  =  l,mdo 

begin 

if  ( tvi  <  tvi>)  then 

Wi  *—  Wi> 

S<-SU{:} 

else 

Wi  *—  Wi 

S'-S'U{t'} 

endif 

end 

(comment:  This  completes  the  selection  of  indices  according  to  Proposition  2.) 
begin 

solve  MP(At,  w)  to  get  D  and  D'\  (comment:  Corollary  3.) 

S  +-  S\D 

s' <- s' \iy 

return  5,  S',  D,  and  D' 
end 
end 

Thus,  as  a  result  of  solving  MP(AT,  w),  one  obtains  the  index  set  E(w)  and  its  parti¬ 
tioning: 

E(w)  =  D  U  D'  U  S  U  S'. 

Indeed,  as  will  be  seen  in  the  next  section,  it  is  useful  to  partition  the  full  sets  of  indices 
as  follows: 

{l,...,m}  =  DU  SU  N  and  =  D'US'UN'. 

The  elements  of  N  and  N'  are,  of  course,  the  indices  of  nonbasic  columns  of  A1.  There  is 
a  pairing  between  the  elements  of  5  and  those  of  N'.  Likewise,  there  is  a  pairing  between 
the  elements  of  S'  and  those  of  N, 
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In  the  example  used  in  Section  3,  these  six  sets  are 


B={2,3},  B’={r,S), 

5  =  {1},  5'  =  0, 

N  =  0,  N’=  {1'}. 

5.  Solving  systems  of  equations. 

Each  iteration  of  the  Bax  Method  entails  the  solution  of  a  box  problem,  each  of  which, 
in  turn,  calls  for  the  solution  of  two  systems  of  equations  stated  in  terms  of  the  current 
minimum-weight  basis.  [See  (12)  and  (13)  in  Zikan  and  Cottle  (1987).]  This  section  will 
discuss  how  to  treat  these  systems  when  the  original  primal  problem  is  of  the  bounded 
variable  type. 

In  general,  let  BT  denote  the  minimum- weight  basis  (A1).^).  The  first  system  of 
equations  to  be  solved  is  of  the  form 

Btc  =  c,  (9) 

where 

•-(:)■ 

The  second  system  of  equations  to  be  solved  is  of  the  form 

Bz  =  z  (10) 

where  z  is  the  solution  to  a  simplified  box  problem  [see  (14)  in  Zikan  and  Cottle  (1987)] 
that  is  defined  in  terms  of  c.  The  solution  z  of  (10)  above  is  the  solution  to  the  box  problem 
corresponding  to  the  current  minimum-weight  basis. 

Up  to  a  rearrangement  of  rows  and  columns,  the  matrix  AT  can  be  written  in  the  form 

f  At.d  At.s  At.n  0  0  0  y 

Id  0  0  Id>  0  0 

0  Is  0  0  IN.  0 

^  0  0  In  0  0  Is>  j 


After  further  rearrangement,  it  follows  that 


bt= 


A\d  0  A\s  0 

Id  Id1  0  0 

0  0  Is  0 

0  0  0  75- 


Note  that  the  submatrix  At.d  of  Br  is  the  minimum-weight  basis  matrix  that  would 
have  been  obtained  if  there  had  not  been  any  explicit  upper  bounds.  Prom  the  (essentially 
block  triangular)  structure  of  2?T,  it  is  evident  that  solving  systems  of  equations  like  those 
indicated  in  (9)  and  (10)  requires  little  more  work  or  storage  than  than  solving  equations 
(12)  and  (13)  in  Zikan  and  Cottle  (1987). 

To  see  this  more  clearly,  consider  a  system  of  the  form  (9).  In  more  detail  the  system 
can  be  written  as 


At.d  0  A\s  0 

Id  Id1  0  0 

0  0  7S  0 

0  0  0  75- 


Solving  (11)  can  be  done  by  sequentially  solving  the  following  set  of  four  smaller  systems: 


Is'Cs'  =  cs<; 

(12) 

7sC5  =  C5; 

(13) 

At.d£d  =  cd  —  At,sCs\ 

(14) 

Id'Cd1  =  cp'  —  Ipcp- 

(15) 

In  block  form,  (10)  becomes 


Id  0 
Id'  0 


the  technique  for  solving  it  is  analogous  to  that  for  (11),  and  will  not  be  spelled  out. 

Remark.  Due  to  the  possible  need  for  scaling,  it  may  happen  that  some  of  the  blocks 
identified  here  as  identity  matrices  are  actually  nonsingular  diagonal  matrices.  It  is  clear 
that  if  these  conditions  obtain,  there  is  very  little  extra  storage  or  computational  effort 
required.  The  only  significant  work  stems  from  equation  (14). 
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The  procedure  presented  here  can  be  modified  to  handle  the  case  where  some,  but  not 
a 11  of  the  variables  are  explicitly  upper  bounded.  The  details  are  easily  worked  out  and 
hence  are  omitted. 

References 

G.B.  Dantzig,  Upper  bounds,  secondary  constraints,  and  block  triangularity  in  linear  pro¬ 
gramming,  Econometrica  23  (1955)  174-183. 

G.B.  Dantzig,  Linear  Programming  and  Extensions  (Princeton  University  Press,  Princeton, 
New  Jersey,  1963. 

K.  Zikan  and  R.W.  Cottle,  The  box  method  for  linear  programming:  Part  I  -  Basic  Theory, 
Technical  Report  SOL  87-6,  Department  of  Operations  Research,  Stanford  University,  June 
1987. 


!i 


» 


REPORT  DOCUMENTATION  PACE 


W 


SOL  87.-9 


«.  T1TUI  (m4 


The  Box  Method  for  Linear  Programming: 

Part  II  —  Treatment  of  Problems  in 

Standard  Form  with  Explicitly  Bounded  Variables 


■  autmoatw 


Karel  Zikan  and  Richard  W.  Cottle 


•.  sanroMiiNa  onoani  cation  nams  an 

Department  of  Operations  Research 
Stanford  University 
Stanford,  CA  94305 


II.  CONTNOLLINO  OMICI  NAM*  ANO  AOOMSS 

Office  of  Naval  Research  -  Dept,  of  the  Navy 
800  N.  Quincy  Street 
Arlington,  VA  22217 


IT.  OtSTNISUTION  (TATKMCNT  (•!  *• 


DO  1  j22*T»  1473  UtTION  Of  I  NOV  M  IS  OMOi.IT* 


MCUMTV  Ck 


MCUWTY  Cl.*— IFtCATtOW  or  THIS  PAMTVfcai 


SOL  87-9:  The  Bax  Method  for  Linear  Programming:  Part  II  —  Treatment  of  Problems  in 
Standard  Form  with  Explicitly  Bounded  Variables,  Karel  Zikan  and  Richard  W. 
Cottle  (July  1987, 11  pp.) 

A  crucial  aspect  of  the  Box  Method  for  linear  programming  is  the  finding  of  a  "minimum-weight 
basis"  "corresponding  to  a  given  interior  feasible  point.  This  subproblem  leads  to  the  formation  of 
the^Box  Problem,"  a  special  linear  program  having  a  closed  form  solution  which  provides  the  search 
direction  at  the  current  iteration.  Finding  a  minimum-weight  basis  is  a  matroidal  (or,  combinatorial) 
optimization  problem  that  can  be  handled  by  a  greedy  algorithm.  This  paper  suggests  a  way  of  efficiently 
solving  the  minimum-weight  basis  problem  in  cases  where  the  (primal,  standard  form)  linear  program 
contains  explicitly  bounded  variables.  It  is  shown  that  the  main  part  of  the  task  requires  almost  no  more 
computational  effort  or  storage  space  than  does  a  problem  of  the  same  size  without  upper  bounded 
variables.  While  this  result  is  believed  to  be  valuable  in  its  own  right,  there  is  additional  benefit  to 
be  gained  in  applications  where  the  finding  of  a  minimum-weight  basis  (for  a  linear  program  without 
explicit  upper  bounds  on  its  variables)  is  done  by  a  special  greedy  algorithm.  Such  is  the  case  with 
minimum-cost  network  flow  problems  which  will  be  discussed  in  Part  III  of  this  series. 


